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ABSTRACT 

We present results of high time resolution photometry of the eclipsing pre-cataclysmic 
variable NN Ser. NN Ser is a white dwarf/M dwarf binary with a very low-mass 
■ secondary star (~ 0.2 M Q ). We observed 13 primary eclipses of NN Ser using the 

high-speed CCD camera ULTRACAM and derived times of mid-eclipse, from fitting 
of light curve models, with uncertainties as low as 0.06 s. The data show that the period 
ly-j ' of the binary is decreasing, with an average rate of P — (9.06±0.06) x 10~ 12 s/s, which 

c"^) • has increased to a rate of P = (2.85 ± 0.15) x 10 -11 s/s over the last 2 years. These 

rates of period change appear difficult to reconcile with any models of orbital period 
change. If the observed period change reflects an angular momentum loss, the average 
loss rate (J = 1.4±0.6 x 10 35 ergs) is consistent with the loss rates (via magnetic stellar 
wind braking) used in standard models of close binary evolution, which were derived 
from observations of much more massive cool stars. Observations of low-mass stars 
such as NN Ser's secondary predict rates of ~100 times lower than we observe. The 
alternatives are either magnetic activit y-driven changes in the quadrupole moment 
of the secondary star (Applegate, ,1992) or a light travel time effect caused by the 



O 
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presence of a third body in a long (~ decades) orbit around the binary. We show that 
Applegate's mechanism fails by an order of magnitude on energetic grounds, but that 
the presence of a third body with mass 0.0043 M < M 3 < 0.18 M Q and orbital period 
30 < P3 < 285 years could account for the observed changes in the timings of NN Ser's 
mid-eclipses. We conclude that we have either observed a genuine angular momentum 
loss for NN Ser, in which case our observations pose serious difficulties for the theory 
of close binary evolution, or we have detected a previously unseen low-mass companion 
to the binary 

Key words: binaries: eclipsing - stars: individual: NN Ser - stars: evolution - stars: 
fundamental parameters - planetary systems 



1 INTRODUCTION 

The evolution of close binary systems is governed by angular 
momentum (AM) loss, driven by a combination of gravita- 
tional radiation (iKraft et alJll962t lFaulknerlll97lft . which 
is dominant for periods P or b < 3h, and magnetic braking 
iVerbunt fc Zwaarl ll981k which dominates for P or i, > 3h. 
Gravitational radiation is relatively well understood, with 
the angular momentum loss rates simply governed by the 
masses and separation of the components of the binary sys- 
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tem. Magnetic braking, however, is a more complicated is- 
sue. The mechanism is driven by the magnetic field and stel- 
lar wind of one of the binary components. As mass is driven 
off in the stellar wind, the ionised particles are forced to 
co-rotate with the field lines out to the Alfven radius. This 
draws angular momentum away from the star, effectively ex- 
erting a braking force to slow its spin. In close binaries, the 
donor star is tidally locked to the primary, so the angular 
momentum loss cannot act to slow the spin period of the 
secondary alone. Instead, the angular momentum is drawn 
from the binary orbit, causing it to shrink and the orbital 
period of the binary to decrease. The rate of angular mo- 
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mentum loss by magnetic braking is governed by the mass, 
radius and angular momentum of the magnetically active 
star, but there is also evidence to suggest that the AM loss 
rate saturates for lo w-mass stars ab ove a certain value of 
angular momentum JSills et al.ll2000t) . This has led to a ma- 
jor revision in magnetic braking loss rates for binaries with 
low-mass secondaries, such as cataclysmic variables. 

Cataclysmic variable stars (CVs) are mass-transferring 
binary systems comprising a white dwarf primary and a low- 
mass main-sequence secondary. The majority have orbital 
periods between 1.3 and 10 hours, and their evolution is 
governed, as with all clos e binaries, by angular momentum 
loss from the system (see IWarnerlll995l for a comprehensive 
review). CVs are useful for testing close binary evolution, 
as any theory is constrained by 2 major features in the dis- 
tribution of CV periods: the period gap and the value of 
the period minimum. The period gap is a dearth of systems 
with periods between 2 and 3 hours. Standard CV theory ex- 
plains this gap by assuming that magnetic braking is cut off 
at P ~ 3 h as the secondary becomes fully convective (there 
is no longer a radiative/convective boundary to anchor the 
magnetic field, so it either dissipates or is rearranged, re- 
sulting in a lowered stellar wind). Mass transfer ceases until 
the system evolves to a period of ~ 2 h and gravitational 
radiation becomes strong enough to drive mass transfer, re- 
populating the period distribution below the gap. 

The value of the minimum period is governed entirely by 
the angular momentum loss rate for short-period systems. 
Under the standard theory, systems below the period gap 
are driven by gravitational radiation alone, which implies 
that the minimum period should be at P m in ~ 1.1 h. In fact 
the observed cut-off is at about P m in — 1.3 h, suggesting 
that gravitational radiation alone is not strong enough to 
reproduce the observed value of Pmin- 

The standard model was developed by extrapolation 
from studies of braking rat e s of solar-type stars in clus- 
ters (iRappaport et al 1 119831: ISpruit fc Ritteil Il983ft . How- 
ever, a recent dramatic increase in t he amount of data 
available for stars in these clusters fsee lstairfferet alj ll997l : 
iKrishnamurthi et alJ Il997l : iReid fc Mahonevl feOOoT for re- 
views) has shown that this extrapolation to lower masses 
appears to be totally wrong - low-mass stars retain more 
of their angular momentum than their higher-mass coun- 
terparts. This means that the new suggested J is anything 
between 10 and 10 4 times smaller than assumed in the ma- 
jority of CV studies. Importantly, there is also no evidence 
for a cut-off in magnetic braking as th e secondary becomes 
fully convective dAndronov et alJ [2003'l , so the new data of- 
fer no explanation for the existence of the period gap. 

We therefore need a way to directly measure the angular 
momentum loss rates of CVs in order to test the standard 
vs. reduced magnetic braking models. One way of doing this 
is to measure mid-eclipse timings of eclipsing binary sys- 
tems to find the period change of the system and calculate 
the angular momentum loss. The period changes are tiny - 
of order 5 x 10 -4 seconds per year - so this is difficult to 
do in CVs due to contamination of the light curve by the 
accretion processes. In order to overcome this problem, we 
have used a non mass-transferring pre-CV, NN Ser. NN Ser 
is ideally suited for this study as, in addition to the uncon- 
taminated light curve, the system also displays deep primary 
eclipses that give rise to very sharp ingress and egress fea- 



Table 1. System parameters of NN Ser. RD : 
star 



: M dwarf secondary 





Catalan et al. 1994 
Wood & Marsh, 1991 


Haetner et al. 2004 


Binary sep. 


0.95 ± 0.025 R 


0.9543 ± O.O233R 


Inclination 


84° < i < 90° 


84.6°± 1.1° 


Mass ratio 


0.18 < q < 0.23 


0.2778 ± 0.0297 


WD mass 


0.57 ± 0.04 M 


0.54 ± 0.05 M 


RD mass 


0.1 < M < 0.14 


0.150 ± 0.008 M 


WD radius 


0.017 < R Q < 0.021 


0.0189 ± 0.0010 R 


RD radius 


0.15 < R < 0.18 


0.174 ± 0.009 R 


WD temp. 


55000K ± 8000K 


57000K ± 3000K 


RD temp. 


2900K ± 150K 


2950K ± 70K 


RD irr. temp. 


5650K < T < 8150K 


7125K ± 200K 


RD spec, type 


M4.7 - M6.1 


M4.75 ± 0.25 


Distance 


356 pc < d < 472 pc 


500 ± 35 pc 



tures. Given the high-time resolution of ULTRACAM, we 
are able to estimate the times of mid-eclipse to an accuracy 
of ~ 0.15 s. 

NN Ser is a white dwarf / M dwarf binary system with 
an extremely low-mass (M ~ 0.15 M ), and therefore fully 
convective s econdary star. The system was first studied 
in detail bv lHaefnerl jl989T) . who identified it as a deeply 
eclipsing (> 4.8 mag) pre-cataclysmic variable with a strong 
reflection effect of ~ 0.6 mag, and an orbital period of 
0.13 days. I Wood fc Marshl <199lft used low-resolution IUE 
spectra to derive the system param eters, which w e re re- 
fined by the radial velocity study of ICatalan et al.1 <ll994t) 
to give the values i n Table Q The most recent study by 
lHaefner et alJ J2004T) combines high-speed photometry from 
the Multi-Channel Multi-Colour Photometer (MCCP) with 
VLT trailed photometry and phase-resolved spectroscopy. 
This allows them to put good constraints on the temper- 
ature of the secondary star. They also attempt to derive 
accurate values for the radii and masses of the system com- 
ponents, but they failed to detect the secondary eclipse for 
NN Ser, which caused them to derive a binary inclination 
of i — 84.6° ± 1.1°. We have detected the secondary eclipse 
m our ULTRACAM data (see Fig.0, and our preliminary 
modelling indicates that the true inclination is i ~ 88° . Full 
results of our modelling will be the subject of a future paper, 
but w e conclude from our initial results that lHaefner et all 
i2004t) have overestimated the radius (and therefore the 
mass) of the secondary star by ~ 15%. Nevertheless, we 
carry out all of our analysis for all val ues of mass and radius 
incl uded in the un c ertain ties given bv lCatalan et all <ll994t) 
and lHaefner et al ] J2004I) . All four studies also give eclipse 
timings (listed in Table[3Jl which we have used to extend our 
baseline for measuring the period change to ~ 15 years. 



2 DATA ACQUISITION 

The data were taken with the ultra-fast, triple-be am CCD 
camera, ULTRACAM fsee lDhillon fc Marshl fcOOll) for a re- 
view). We used the camera in conjunction with the 4.2 m 
William Herschel telescope at the INC to observe NN Ser 
simultaneously in the Sloan u', g' and either r', i' or z' bands. 
We set a time resolution of ~ 2.06 s to ensure a high signal- 
to-noise ratio in all wavebands. The observations were taken 
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Table 2. 

NN Ser 



Observation log for ULTRACAM observations of 



Date 


Filters 


No. 


Conditions 






eclipses 




17/05/2002 


u'gV 


2 


good, seeing ~1.2" 


18/05/2002 


u'g'r' 


1 


variable, seeing 1.2" - 2.4" 


19/05/2002 


u'g'r' 


1 


fair, seeing ~2" 


20/05/2002 


u'g'r' 


1 


fair, seeing ~2" 


19/05/2003 


u'g'z' 


1 


variable, seeing 1.5" - 3" 


21/05/2003 


u'g'i' 


1 


excellent, seeing ~1" 


22/05/2003 


u'g'i' 


1 


excellent, seeing <1" 


24/05/2003 


u'g'i' 


1 


good, seeing ~ 1.2" 


03/05/2004 


u'g'i' 


3 


variable, seeing 1.2" - 3.2" 


04/05/2004 


u'g'i' 


1 


variable, seeing 1.2" - 3" 



over a period of 2 years, in May 2002 - 2004, and during 
those runs we observed 13 primary eclipses of the system. 
We were also able to observe a number of secondary eclipses 
as the white dwarf transited the secondary star (see Fig. [Q. 
The pixel size for the 3 ULTRACAM CCDs is 13 /im, with a 
scale of 0.3"/pix. Readout noise is 3.10 - 3.40 e, depending 
on the CCD, while the gain is 1.13 - 1.20 e/ADU.Each UL- 
TRACAM frame is time-stamped to a relative accuracy of 
better than 50 /xs and an absolute accuracy of better than 
1 millisecond using a dedicated GPS system. A full list of 
observations and observing conditions is given in Table [2] 
The weather in May 2004 was particularly variable, leading 
to larger errors in our measured times for that epoch. 

The data were reduced using the ULTRACAM pipeline 
software written by TRM. Differential photometry was per- 
formed on the target, with respect to a nearby, bright, non- 
variable comparison star. 



3 ANALYSIS & RESULTS 

All MJD times were corrected to Barycentric Dynamical 
Time (TDB), then additionally corrected for light travel 
time to the solar system barycentre. All times are there- 
fore listed in MJD(BTDB). In order to measure accurate 
eclipse times, we needed a model of the eclipse of the white 
dwarf, which we calculated as follows. We defined the two 
stars by their radii relative to the separation of the binary. 
Since we allowed for tidal deformation of the M dwarf (but 
not the white dwarf), the radius of the M dwarf was mea- 
sured from its centre of mass towards the white dwarf. Apart 
from the relative radii, we also require the binary mass ra- 
tio and inclination, stellar effective temperatures and linear 
limb darkening coefficients to define our model binary. The 
two stars were divided into many small elements. The tem- 
peratures of the elements covering the M dwarf were set, 
accounting for incident flux from the white dwarf by adding 
fluxes so that 

tT' 2 = + firr, 

where a is the Stefan-Boltzmann constant and Fin is the 
flux incident on the secondary accounting for the projection 
effects and the distance from the white dwarf. The surface 
brightness of each element was then set assuming black-body 
spectra, and given the effective wavelength of the filter in 
question. Once the surface brightnesses were set, the model 




1.95 2 2.05 

MJD(BDTB) - 53410 (days) 



Figure 1. Differential light curves for NN Ser, taken simultane- 
ously in the u', g' and r' Sloan niters from top to bottom respec- 
tively. The light curves are binned by a factor of 43. The hump in 
the light curve is caused by the reprocessing of light from the WD 
by the cool secondary star. A shallow secondary eclipse can be 
detected at the top of the reflection hump in the r' and g' bands. 
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Figure 2. Differential light curve for NN Ser in g' with light curve 
model overplotted. 



light-curves were computed by summing over all elements, 
testing for which were in view and not eclipsed and account- 
ing for their projected areas. The eclipse by the M dwarf 
was computed, allowing once again for tidal distortion. Our 
assumption of black-body spectra for the two stars is physi- 
cally unrealistic, but for the eclipse times of this paper, the 
key element is to have a model that can match the shape 
of the primary eclipse, which ours does well (Fig. The 
timings and associated errors for the mid-eclipses in all of 
the wavebands are given in Table El The errors on our mid- 
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Table 3. Measured times of mid- eclipse for eac h of the 13 observed primary eclipses of NN Ser. Cycle numbers are counted from first 
measured eclipse in the literature jHaefner 1989|. Times were measured for all 3 wavebands simultaneously observed with ULTRACAM 
and are given in MJD(BTDB), i.e. MJD shifted to the solar system barycentre and corrected for light travel time. The red filter varied 
between nights, so the filter used is listed in the final column. The poor observing conditions during eclipse cycle 44474 led to the loss 
of data in the u' band. 



Cycle 


u' eclipse 


uncertainty 


g' eclipse 


uncertainty 


r'/i'/z' eclipse 


uncertainty 


red 


Number 


time 


1 cr 


time 


1 cr 


time 


1 cr 


filter 


38960 


52411.9470588 


0.0000020 


52411.9470564 


0.0000005 


52411.9470577 


0.0000010 


r' 


38961 


52411.0771385 


0.0000016 


52412.0771385 


0.0000005 


52412.0771383 


0.0000010 


r' 


38968 


52412.9876761 


0.0000030 


52412.9876977 


0.0000008 


52412.9876721 


0.0000013 


r' 


38976 


52414.0283427 


0.0000030 


52414.0283394 


0.0000006 


52414.0283379 


0.0000016 


r' 


38984 


52415.0689716 


0.0000025 


52415.0689810 


0.0000007 


52415.0689795 


0.0000016 


r' 


41782 


52779.0331646 


0.0000021 


52779.0331696 


0.0000010 


52779.0331362 


0.0000100 


z' 


41798 


52781.1144524 


0.0000015 


52781.1144513 


0.0000006 


52781.1144567 


0.0000014 




41806 


52782.1550904 


0.0000021 


52782.1550929 


0.0000006 


52782.1550948 


0.0000011 




41820 


52783.9762155 


0.0000022 


52783.9762151 


0.0000007 


52783.9762110 


0.0000020 




44472 


53128.9486787 


0.0000070 


53128.9486778 


0.0000040 


53128.9486611 


0.0000800 




44473 


53129.0787555 


0.0000027 


53129.0787597 


0.0000022 


53129.0787487 


0.0000050 




44474 


no data 


n/a 


53129.2088356 


0.0000020 


53129.2088355 


0.0000027 




44480 


53129.9893197 


0.0000050 


53129.9893229 


0.0000025 


53129.9893148 


0.0000040 


i' 



Table 4. Previous eclipse times of NN Ser (1) Haefner (1989); 
(2) Wood & Marsh (1991); (3) Pigulski & Michalska (2002); (4) 
Haefner et al. (2004). 



Time of mid-eclipse 


Cycle 


Reference 


MJD(BTDB) 


Number 




47344.025(5) 





1 


47703.045744(2) 


2760 


4 


47703.175833(6) 


2761 


4 


47704.216460(3) 


2769 


4 


47705.127023(3) 


2776 


4 


47705.257115(7) 


2777 


4 


47712.28158(15) 


2831 


2 


47713.32223(15) 


2839 


2 


48301.41420(15) 


7360 


2 


51006.0405(2) 


28152 


4 


51340.2159(2) 


30721 


4 


51666.9779(4) 


33233 


3 



eclipse timings are typically ~ 0.15 s but as low as 0.06 s 
when conditions are good. 

The times of mid-eclipse in the g' band were plotted 
against cycle number. We found that all 13 of the ULTRA- 
CAM points except for one (cycle 38968) fell on a straight 
line, and that the one discrepant point showed a time shift of 
2.06 s - exactly the same timing as one exposure. We noted 
from the logs that we had GPS problems during this expo- 
sure run, therefore concluded that the GPS timestamp had 
slipped by one exposure for that point, and corrected it by 
2.06 s to bring it in line wit h the other points. Old eclipse 
timings from the literature jHaefnerlll989t IWood fc Marshl 
Il99ll: iPigulski fc Michalskal 120021 lHaefner et al.ll2004l ; Ta- 
ble0J were then added to the plot. The residuals after sub- 
tracting a straight-line fit can be seen in Fig[3] We derive 
a best-fit linear ephemeris from all the available data for 
NN Ser as 

MJD(BTDB) = 47344.0246049(14) + 0.130080144430(36)E, 



where the quoted uncertainties are the 1 sigma uncer- 
tainties in the fit. We derive a best-fit quadratic ephemeris 
as: 

MJD(BTDB) = 47344.0244738(16) + 0.13008017141(17)E 

-5.891(36) x 10~ 13 E 2 . 

The data, with eclipse times t_B , were fit with a parabola 
of form: 

ts = To + AE + BE 2 . (1) 
The rate of period decrease can then then be found using 

(2) 



We found that the rate of period change over the 15 
years of observations is increasing, so we fit all of the data 
to find an average rate of period change, and just the UL- 
TRACAM data to find the current rate. 
The angular momentum of the system as a whole is given 
by 

J =(w) M1M2, (3) 

where Ml, M2 and M are the primary, secondary and total 
masses respectively. Combining this with Kepler's third law, 



4tt 2 



GM 

,3 ' 



(4) 



P 2 a" 

we find that, for a detached system (where Ml, M2 and M 
are constant), 



(•>) 



J _ 2B_ 
J ~ 3P2' 

For NN Ser, our measured value for the average rate of pe- 
riod change is 

Pav = 9.06 x 10~ 12 ± 0.06 x 10" 12 s/s 
and for the current rate of period change 
Pcur = 2.85 x 10~ n ± 0.15 x 10~ n s/s. 
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Taki ng 0.1M W < M 2 < 0.14 M w and 0.15 R ^ R 2 ^ 
0.18 iCatalan et al.lll994T) . these correspond to angular 
momentum loss rates of 



0.! 



x 10 35 < Jav < 2.09 x 10 35 ergs. 



and 



2.52 x 10 35 < j cur < 6.87 x 10 35 ergs. 

where the relatively large allowed range is caused by the 
uncertainties in the system parameters, where we have as- 
sumed that the system parameters are independent of each 
other. Obviously this will overestimate the size of the uncer- 
tainty in our measured angular momentum loss rate for any 
one value of secondary mass. When carrying out the analysis 
of period loss mechanisms in Section 4, we have used a more 
realistic approach, relating the secondary radius to its mass, 
us ing the M-R relation f or sec ondaries in binary stars given 
in lGorda fc Svechnikovl dl998l) and calculating the resultant 
separation of the binary. This then gives a range of values 
for the angular momentum change that are specific to each 
value of secondary mass. 



DISCUSSION 
CHANGE 



MECHANISMS FOR PERIOD 



Period changes in binary systems are generally due to one 
of three mechanisms: 

(i) Applegate's (1992) mechanism, where period changes 
are caused by coupling between the binary period and 
changes in the shape of the secondary star; 

(ii) the presence of a third body in a long orbit around the 
binary. This affects the light travel-time, which can be mis- 
interpreted as a change in the binary period. For example, as 
the binary moves towards the observer, the eclipses are seen 
to occur more frequently than when the binary is moving 
away. 

(iii) a genuine angular momentum loss from the system; 

We show below that the most common cause of mea- 
surable change in binary periods - Applegate's mechanism 
- cannot work for NN Ser: the luminosity of the secondary 
star is too low to provide the necessary energy. We also 
discuss the other two mechanisms in detail, along with the 
ramifications for binary evolution. 



4.1 Applegate's mechanism 

lApplegatd" dl992T) proposed that orbital period modulations 
observed in many binary stars could be induced by the 
gravitational coupling of the binary orbit to variations in 
the shape of the magnetically active secondary star. The 
shape changes are reflected in a change of quadrupole mo- 
ment which leads to the change in period; no loss of an- 
gular momentum from the system is necessary. The shape 
changes are presumed to be driven by solar -like magnetic cy- 
cles. T o avoid an excessi ve energy budget iMarsh fc Pringld 
ll990T) . lApplegatd" il992l) proposed that the shape changes 
were driven by the re-distribution of angular momentum 
within the secondary star. He showed that the energy re- 
quired was well within the capabilities of 4 out of 5 stars 



that he considered, and was not far off the mark for the 
fifth, RS CVn. 

One of the reasons we chose NN Ser was that it has 
a particularly low mass - and therefore low luminosity - 
secondary star, which should be i ncapa ble of driving large 
period changes under lApplegatd <ll992T) 's model. We now 
consider this in detail. 

The observational fact we have to explain is the to- 
tal period change, which in NN Ser amounts to AP = 
(-4.26 ±0.03) x 10~ 3 s over the 15 years from the MCCP to 
the ULTRACAM epoch. A period change AP corresponds 
to a change in quadrupole moment AQ, where 

AP (R\ 2 AQ 

where M and R are the mass and radius of the secondary 
star and a is the orbita l separation iApplegate fc Pattersonl 
Il987h . lApplegatd 11992T) calculated the change in quadrupole 
moment by considering the transfer of angular momen- 
tum from the inside of the star into a thin outer shell. 
This increases the oblateness of the shell and therefore 
its quadru pole moment, at the expense of some energy. 
lApplegatd" used M s ~ 0.1 Mq for the mass of the shell. 
This is immediately a problem in the case of NN Ser be- 
cause the mass of the secondary star is only O.15M0, and 
so 0.1 M(?) is not in any sense a "shell". We therefore gen- 
eralised ^^2^£^3' s work as follows. We split the star into 
an inner "core", denoted by subscript 1 and an outer shell 
denoted by subscript 2. Angular momentum is transferred 
from the core to the shell leading to changes in their angular 
frequencies of AQi and Afi 2 such that 



ZiAfii + hAn 2 = 0, 



(7) 



where I\ and I2 are the moments of inertia and given by 
integrals over radius of a series of thin shells of the form 



I = ~ / p-'d.W. 



(8) 



The change in angular frequencies changes the oblateness 
and therefore quadrupole moments of the core and shell lead- 
ing to an overall change in quadrupole moment of 

AQ = Qi [2fii Afii + (Afii) 2 ] +Q 2 [2fi 2 Afi 2 + (Afi 2 ) 2 ] , (9) 

where the Q' coefficients are given by integrals over shells of 
the form 



Q' 



1 J R b dM 
9 J GM(R) ' 



(10) 



where M(R) is the mass inside radius R. These equations 
follow from Eq. 25 of lApplegatd" dl992l) . 

For a given period change, mass and radius Eq.|S] gives 
the change in quadrupole moment. We then use Eqs[7|and|5| 
to solve for AS7 2 , and therefore for the angular momentum 
transferred, AJ = I 2 Af2 2 . T his then leads to the energy 
change from lApplegatd <1992l) 's Eq. 28: 

1/1 1 \ 1 

(AJ) 2 , (11) 



ap = n dr Aj + 

2 Vii 12 

where Qdr = f2 2 — fii is the initial differential rotation. 

In order to solve the equations above, one needs first 
to know the run of density with radius. We calculated this 
from the Lane-Emden equation for an n — 1.5 polytrope as 
an approximation to the fully convective secondary star. 
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Figure 3. The upper plot is an O-C diagram showing the period change in NN Ser. A linear fit has been subtracted from the data. 
The solid line is a fit through all the data (average rate of period change), while the dashed line is a fit through the ULTRACAM data 
only (current rate of period change). The lower three plots are (from uppermost): residuals after a fit through all the data is subtracted, 
showing all the points; residuals after the fit through all the data is subtracted, zoomed in on the ULTRACAM points; residuals after 
the fit through the ULTRACAM data is subtracted, zoomed in on those points. 



In Fig. 2] we show the results of these calculations as a 
function of the shell ma ss for M = 0.15 M , R = 0.174 R 
and AP = —0.00426 s. lApplegatel used a value of fW = 
AQ2, on the basis that one would expect the initial differ- 
ential rotation to be of a similar order of magnitude as the 
changes. In order to arrive at a minimum energy, we have 
assumed that Qd r = 0; had we used instead fl^r = Af^, 
the energies would be increased by a factor ~ 2. Even 
without this factor, the figure shows that we need at least 
~ 4 x 10 40 ergs to drive the observed period change. The 
luminosity of the secondary star is L2 = iirR* ''cT^f f, which, 
for NN Ser's secondary with 2880 K < T c ff < 3020 K, gives 
the energy available over the 15 years of observations as 

2.5 x 10 39 ergs < E 2 «S 4.5 x 10 39 ergs. 

This range, which is indicated in the lower-right of Fig. 0] 
fails by a factor of about 10 to match the energy we cal- 
culate, although it does just tally with a calculation based 



upon lApplegatel ' s original calculations (dotted line). A plot 
of the ratio of our value of minimum energy required to drive 
Applegate's mechanism over the energy available in 15 years 
vs the radius of the secondary star is shown in Fig|3 It can 
be seen that the ratio of energy required over energy avail- 
able is well above 1 for all values of secondary star radius 
(a nd therefore mass), i. e. for the system paramet ers derived 
bv lCatalan et alJ Jl994l) and lHaefner et alJ (I2004D . NN Ser's 
secondary star is not capable of generating enough energy 
to drive Applegate's mechanism. 

Our increased estima te is a resu l t of t he differences in 
our approach compared to lApplegatel 1^992) 's. First we gen- 
eralise his thin shell approximation by integration over finite 
ranges of radii. Second, in Eq. [^J there is a positive contri- 
bution from the shell, but also a negative contribution from 
the core which lApplegatel did not include and which bal- 
ances the shell to a considerable extent. The dashed line 
in Fig. 21 shows the effect of ignoring the core's contribu- 
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Figure 4. A plot of the energy requir ed to effect the period 
change in NN Ser using lApplegatd ll992l) 's mechanis m as a func - 
tion of assumed shell mass. The dotted line shows lAppIegatd 's 
original calculation. The solid line shows our calculation, inte- 
grating over shells and allowing for the quadrupole moment of 
the inner core of the star. The dashed line shows the result if 
we ignore the quadrupole moment of the core. The error bar in 
the lower-right shows the energy available to effect the change in 
NN Ser. 



tion to the energy in our calculation. This line converges 
towards lApplegatef s (dotted line) at small shell masses. At 
large shell masses the effect of the var iation in de nsity with 
radius and finite shell thickness, which lApplegatel did not in- 
clude, are important and explain the remaining difference. 
Ignoring the core cuts the energy requirement by about a 
factor of 4, suggesting that, given the squared dependence of 
energy on A J and hence AQ, the core balances about 50% of 
the quadrupole increase from the shell. lApplegatel il992T) al- 
ready recognised that his approximations must break down 
when the shell becomes a significant fraction of the star's 
mass; the counter-balancing effect of the core's quadrupole 
moment has not been pointed out before as far as we are 
aware. 

While we have shown that the intrinsic luminosity of 
NN S er's secondary star is too low to drive lApplegatel 
( 1992) 's mechanism, we note that the secondary star is heav- 
ily irradiated by the hot, white dwarf primary, as shown by 
the much higher temperature on the side facing the white 
dwarf. In order to estimate how much of an effect this will 
have on the secondary star, we compared the flux from the 
primary that is intercepted by the secondary star with the 
secondary star's intrinsic luminosity. We find that the inter- 
cepted flux from the white dwarf is ~ 13 times the intrin- 
sic luminosity of the secondary, and therefore, if more than 
70% of this were to be absorbed by the secondar y star, the 
extra energy provided could be enough to drive lApplegatel 
( 1992) 's mechanism. However, there are a number of reasons 
that this should not be the case. Firstly, very little of the 
energy absorbed at the stellar surface flows inward, since the 
opacity in this region becomes very high, the stellar surface 
quickl y heats up, and so the heatin g luminosity is re-radiated 
(Har paz fc RappaportHl99ll .ll995"). Instead, the main effect 
of the irradiating flux is to block the outflow of the radia- 
tion produced in the secondary star's interior. The star will 
undergo a small expansion as some of the blocked energy 
is stored as internal or gravitational energy, but unlike a 
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Figure 5. Ratio of energy required for Applegate's mechanism 
over energy generated by the secondary star, vs secondary radius. 
The dashed lines show the uncertainty in energy ratio due to the 
uncertainties in the observed period change and the temperature 
of the secondary. The dotted line is at a ratio of 1, i.e. Applegate's 
mechanism for period change is only possible below this line. 



star undergoing isotropic irradiation, an anisotropically ir- 
radiated star diverts the energy flow in the upper layers 
of the convection zone to the unirradiate d parts of its sur- 
face, efficiently cooling t he secondary star llVaz fc Nordlandl 
Il985t iRitter et ai]l2000l) . In this case the energy flow is de- 
coupled from the mechanical and thermal structure, which 
can still be considered as spherically symmetrical, hence the 
structure of the secondary star below the convection zone is 
virtually unaffected by the irradiating flux. The irradiating 
flux does not penetrate deeply enough into the atmosphere 
to affect the deep bounda ry layers that must be deformed 
to drive IXpplegatd il992T) 's mechanism. Finally, if the sec- 
ondary star were absorbing ~ 9 times its intrinsic energy 
from the primary star, we should see some variability in the 
light curves over the 3 years of observations with ULTRA- 
CAM. We see no evidence for this as the light curves are 
consistent to within 3.5% over that time. 

W e conclude that in the case of NN Ser, lApplegatel 
( 1992) 's quadrupolar distortion mechanism falls short of be- 
ing able to match the observed period change, although we 
note that there may yet be w ays to affect the q uadrupole 
distortion at less energy cost llLanza et alJll993) . We now 
look at alternative mechanisms of producing NN Ser's pe- 
riod change. 



4.2 Third body 

Apparent changes in the orbital periods of binary stars 
have often been attributed to the light-travel time variation 
caused by third bodies, although further observation usually 
reveals that this cannot be the case. However with the rela- 
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tively limited coverage to date, this is at least a possibility 
for NN Ser which we investigate in this section. 

Changes in eclipse timings of binary stars do not nec- 
essarily indicate a genuine change in the binary period. A 
third body in a long orbit around the binary can cause small 
but significant changes in the light-travel time from the bi- 
nary system, which manifest themselves as strictly sinusoidal 
changes in the timings of mid-eclipse. We are able to put 
constraints on the mass and period of any third body which 
could cause the observed period change in NN Ser by fitting 
all possible sine waves to a plot of mid-eclipse timings vs 
cycle number. A function of form 

T — To + P orb E + A 3 sin ( ^< E ~ E ^ 

was fitted to the plot for values of P3 between 2 and 500000 
days, where P or b was kept fixed at the orbital period of the 
binary, P3 is the modulation period of the period change, 
A3 is the amplitude of the period modulation, and E — E3 is 
a measure of the phase of the zero point of the modulation 
with respect to the zero point of the binary period. As we 
are interested in the minimum possible mass, we assumed 
that the inclination of the orbital plane of the third body is 
aligned with the line of sight, i.e. sin i = 1. This gave us the 
values of A3 for all possible modulation periods between 2 
and 500000 days. From this, we were able to use Kepler's 
law and the observed luminosity of NN Ser to find the range 
of allowable masses of the third body which could cause the 
observed period change in NN Ser. The minimum possible 
mass comes from the fact that we have not seen a reversal 
in the period change of NN Ser. The minimum value of P3 
is therefore ~ 30 years, which corresponds to a minimum 
mass of M 3 = 0.0043 M Q . 

The maximum value of M3 comes from the luminosity 
of the binary system in eclipse. The luminosity of the third 
body must be equal to or less than the observed mid-eclipse 
luminosity. This means that it must have a mass equal to 
or less than that derived for the secondary star. If the max- 
imum mass is O.I8M0 then the maximum orbital period for 
any third body is P3 = 1.04 x 10 5 days ~ 285 years. 

We therefore find that, on the basis of the current data 
at least, a low-mass companion to the binary system could 
cause the observed changes in mid-eclipse timings that we 
observe in NN Ser, and that the long periods suggested by 
our data would be able to accommodate NN Ser's primary 
even before its evolution to a white dwarf. Our results also 
indicate that measuring eclipse timings of binary systems is 
potentially a very sensitive method of detecting extra-solar 
planets in long-period orbits. We suspect however, that as 
in other instances, further observations will rule out a third 
body. 



4.3 Comparison with angular momentum loss 
models 

The period decrease we have measured in NN Ser may also 
be explained by angular momentum loss from the binary 
system. Angular momentum loss in CVs and pre-CVs is gov- 
erned by two mechanisms - gravitational radiation and mag- 
netic braking. The rates of AM loss caused by both mech- 
anisms must be added together to find the total AM loss 



for the system. We compare the inferred angular momen- 
tum loss rate for NN Ser (corresponding to the rate of pe- 
riod decrease) to both the v alues predicted by the s tandard 
CV magnetic braking rate iRappaport et al.Ml983T) . based 
on extrapolation from studies of braking rates of solar-type 
stars in clusters, and to the reduced magnetic braking rate 
JSills et al . 2000), based on more recent data, for which the 
angular momentum saturates at lower masses. Under the 
standard model, the angula r momentum, J, decreases as 
J ex —lo iSkumanichlll972l) . where lo is the angular ve- 
locity of the star. However, the reduced braking model sug- 
gests that the angular momentum loss is best modelled as 
J oc —lo' 3, for lo < uicrit and J oc — lo for lo > Lo cr u, where the 
threshold rate, Lo cr i t is much lower than the rotation rates 
of CVs. This means that the new suggested J is anything 
between 10 and 10 4 times smaller than assumed in the ma- 
jority of CV studies. 

If this is correct, we require a large-scale revision of 
CV evolution, possibly with systems staying at an approx- 
imately fixed period throughout their lifetime rather than 
migrating from long to short periods. However, such a model 
has significant problems when compared to observations, 
particularly as the mass transfer rate should be much lower 
than seen in the high accretion rate group of CVs known as 
novalike variables. 

Both models w ere applied to CV studies by 
lAndronov et alJ J20031) . hereafter APS03. 

4-3. 1 Gravitational radiation 

We use the same expression for AM loss due to gravitational 
radiation as used in APS03, although this was misquoted in 
their paper. The correct expression is given by: 

/dJ\ = _32G^ a _ 7/2M 2 M 2 Ml/2 (12) 

V dt J grav 5 C 

where Mi, M2 and M are the white dwarf mass, secondary 
mass and total mass respectively, and a is the binary sep- 
aration given by Newton's form of Kepler's third law a = 
(GM/w 2 ) 1 ' 3 . For NN Ser this gives a range of values of 5.75 
xlO 32 ergs < jgrav < 1-74 xlO 33 ergs, over 100 times smaller 
than required to drive our measured value of P for NN Ser. 

4-3.2 Standard magnetic braking model 

The standard model for magnetic braking in CVs is based 
upon studies of the solar wind and a nd the rotation periods 
of solar- type sta r s in open cluster s JWejoexJ^IOavis] ^67|; 
ISkumanichlll972t iMestel fc Sp7uiilll987l) . I Rappaport et alJ 
l^83TT developed an empirical prescription that is still com- 
monly used in CV studies. This relationship is given by 

(^-) « -3.8 x 10~ 30 M Q R%m 2 rluj 3 ergs, (13) 
\ at J m b 

where ^ 7 ^ 4 is a dimensionless parameter and lo is the 
angular frequency of rotation of secondary star (= binary 
period for CVs) in rad s _1 . We applied this to NN Ser to 
find the predicted standard angular momentum loss rate for 
this pre-CV. The results can be seen in Fig.|S| APS03 cut off 
the standard magnetic braking model at a secondary mass of 
0.3 Mq to satisfy the standard CV theory. This states that 
as the secondary becomes convective, the magnetic field is 
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no longer locked to the stellar core and so either dissipates or 
is rearranged, cutting off the magnetic braking mechanism. 
APS03 suggested that there is no evidence for this cut-off, 
so we have not applied it here. We find that by ignoring 
the magnetic braking cut-off, this model can explain the 
observed loss rates seen in NN Ser. 



4-3.3 Reduced magnetic braking model 

The more recently proposed model for angular momen- 
tum loss due to magnetic braking was applied to CVs 
in AP S03. Studies of th e rotational periods of low-mass 
stars (iQueloz et all Il99fil ICollier-Cameron fc Jiankeill994t 

EKcppe ns et al .1 1199.4 iKrishnamurthi et al.lll997t ISills et alJ 
2000) all showed that the standard model overestimates an- 
gular momentum loss rates for periods below 2.5 - 5 days 
and that a modification of the standard model was required 
for those high rotation rates. APS03 modelled the modified 
angular momentum loss rates u sing a pres c riptio n with the 
same functional form as that of lSills et alJ <|2000j). given by 



for Ld < LUcrit 
for UJ > LUcrit 



(14) 



where uj cr it is the critical angular frequency at which the 
angular momentum loss rate enters the saturated regime. 
The constant K w = 2.7 xlO 47 g cm s _1 is calibrate d to give 
the k nown solar rotation rate at the age of the Sun llKawaler] 
Il988l) . The values of io cr it were calc ulated from the values 
of LOaritQ given in ISills et all i2000T) using the relationship 
between u cr it and convective turnover time, r given by 



T0 

LUcrit — LOcritQ • 



(15) 



The values of r were taken from iKim fc Demarauel l(l£ 
assuming an age of 0.2 Gyr. 

Again, the prescription was applied to NN Ser. Results 
are shown in Fig. |S] Our plot differs significantly from the 
original plot in APS03 due to their mis-calculation of the 
angular momentum loss due to gravitational radiation. By 
applying the correct gravitational radiation loss rate, the 
total angular momentum loss rates predicted by APS03 are 
within ~ 1 order of magnitude of the standard magnetic 
braking model rather than ~ 2 orders of magnitude lower, 
as they originally claimed. However, this is still too low to 
explain the loss rates seen in NN Ser. 



5 DISCUSSION 

We have found that only two mechanisms can explain the 
observed period change in NN Ser - either a genuine angular 
momentum loss from the system or an unseen third body in 
orbit around the binary. In the case of an angular momen- 
tum loss, our observations show that the s ystem is losing an- 
gular momentum at the rate predicted bv lRappaport et alJ 
( 1983), but only if we assume that magnetic braking is not 
cut off as the secondary mass reaches 0.3 M Q . APS03 pointed 
out that an increase in the angular momentum loss rate at 
low periods can solve a major problem regarding the the- 
oretical vs observed values of the period minimum. If CV 
evolutio n at P < 2h wer e driven solely by gravitaional ra- 
diation, IPattersonl |l984) found that the period minimum 
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Figure 6. Plot of the braking rates predicted by gravitational ra- 
diation and the standard and reduced braking models for NN Ser. 
The different plots for the standard model are for different values 
of 7 = 0,1,2,3,4, from the top down. The shaded region shows our 
measured value of the braking rate of NN Ser. 



should be at 1.1 hr rather than the observed value of 1.3 hr. 
He also noted that the angular momentum loss rates would 
be very low for these short-period systems, implying a low 
mass accretion rate and therefore a high population of CVs 
at the minimum period, a prediction that is contradicted by 
observation. By adding the extra angular momentum loss 
rate from magnetic braking, the cut-off is shifted to longer 
periods. However, our value of the magnetic braking rate 
causes the opposite problem. We find that at short peri- 
ods, the magnetic braking rate is almost 100 times the rate 
of angular mom entum loss due to gravitational radiation. 
iPaczvnskH jl98lT) showed that, for constant j, 



_J_ 

J or 



(see also |Pattersorj [l984 ) . bringing the minimum period up 
to a value of 331 mins ~ 5.5 hours, which is clearly not 
correct. 

We also have the continuing problem of how to explain 
the presence of the period gap - a dearth of systems with 
periods 2 h P ^ 3 h. If magnetic braking does not shut 
off as the secondary becomes fully convective then there is 
no reason for systems to cease mass-transfer between those 
periods. APS03 suggested that instead of a migration of CVs 
from long- to short-period, the systems above and below the 
period gap may belong to two different populations, with 
no migration between the two. However, this is more likely 
for their longer-timescale angular momentum loss, as their 
model depends upon the presence of an evolved secondary 
star, than for our measured magnetic braking rate. 
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6 CONCLUSIONS 

We find that there are two possible explanations for the ob- 
served period change in the pre-CV NN Ser over the last 15 
years. If the change is due to a genuine angular momentum 
loss from the system then the rate corresponds to an angu- 
lar momentum loss that agrees most close ly with the stan- 
dard magnetic braking rate proposed by IRappaport et all 
l|l983j), and that the reduced magnetic braking rate of 
lAndronov et al.l <l2003f) underestimates the measured rate 
by ~ 2 orders of magnitude. We find no evidence for a 
cut-off in magnetic braking as the secondary mass drops 
below M = 0.3 Mq. If the period change is instead due 
to a third body, we place constraints on such a body of 
0.0043 M < M 3 < 0.18 M Q and 30 years < P 3 < 285 years. 

As a by-product of this investigati on, we hav e found 
that the energy requirements of the lAppleeatel (Il992l) 
quadrupolar distortion mechanism are significantly in- 
creased once one accounts for the role of the inner part of 
the star in counter-balancing the outer shell. 
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